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Abstract 

In this paper we continue the investigation of the Maxwell-Landau-Lifschitz and Maxwell-Bloch equations. In 
particular we extend some previous results about the Cauchy problem and the quasi-stationary limit to the case 
where the magnetic permeability and the electric permittivity are variable. 
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1 Introduction 



The models. This paper deals with two physical models which describe the propagation of electromagnetic waves, 
that is of the magnetic field H and of the electric field E, in some special medium which occupies an open subset of 
M 3 , with magnetic permeability p and electric permittivity e. In both cases we denote by / the extension of a function 
/ by outside the set f2. The time variable is t ^ 0, and the space variable is i £ R 3 . 

The first model refers to Maxwell-Landau-Lifschitz equations (see [10] and [26] for Physics references). The 
magnetic field H and the electric field E satisfy the Maxwell equations in R 3 : 

{pd t H + cur\E 
edtE — curl if 
div p(H + M) 
div sE = 0, 

where M stands for the magnetic moment in the ferromagnet ft and takes values in the unit sphere of Mr . It is solution 
to the Landau-Lifschitz equation: 

d t M = jM AH T - aM A (M A H T ) for x £ 0, (2) 

where 7 =/= is the gyromagnetic constant, and a > is some damping coefficient. Neglecting the exchange phe- 
nomenon, the total magnetic field Ht is the sum 

H T = H + H a (M) + H cxt , (3) 

where the anisotropy field writes H a (M) = Vm$(M)i for some convex function $, and H cyLt is some applied (exterior) 
magnetic field. 

The second model refers to Maxwell-Bloch equations (see for example [7], [8], [16], [33], [36], [39]). In this 
setting Q, denotes some quantum medium with N <G N energy levels described by a Hermitian, non-negative, N x N 
density matrix p. Assuming the usual dipolar approximation, these quantum states change under the action of an 
electric field E by the quantum Liouville-Von Neumann (or Bloch) equation: 

id tP = [A-E-T,p] + iQ(p). (4) 

The N x TV Hermitian symmetric matrix A, with entries in C, represents the (electromagnetic field-) free Hamiltonian 
of the medium. The dipole moment operator T is a N x N Hermitian matrix, with entries in C 3 , and depends on the 
material considered. The (linear) relaxation term Q(p) takes dissipative effects into account (see [5], [5], [3D]). The 
polarization P of the matter is given by the constitutive law P = Tr(Tp) which influences back the electric field E. 
Again, the electromagnetic field satisfies the Maxwell equations in R 3 : 

' pd t H + cm\E = 0, 

ed t E-cmlH = -d t P, 

— (5) 
div(e£ + P) = 0, 

^ div pH = 0. 

Cauchy problems. We first address the questions of global existence, uniqueness and stability for the Cauchy 
problem associated with these equations. The physically relevant solutions have finite energy: they satisfy the usual 
(I/ 2 ) energy estimates. Mathematically, this regularity leads to weak solutions and is usually not enough to ensure the 
desired uniqueness and stability properties (requiring for these hyperbolic semilinear systems in space dimension 3, in 
the general theory, H s Sobolev regularity with s > 3/2). 

However, in the case of the Maxwell-Landau-Lifschitz system, Joly, Metivier and Rauch [23] noticed that specific 
(algebraic) properties of the nonlinearities, as well as (geometric) properties of the differential operator involved, 
allowed to show the existence of global finite energy solutions (essentially, using compensated compactness arguments) 



-pd t M, 



0. 

= 0, 



(1) 
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enjoying stability properties. Furthermore, only a small amount of regularity (curl H and curli? in L 2 ) ensures 
uniqueness. This is achieved using dispersive properties of the system; namely, a limit Strichartz estimate controlling 
the L 2 L<^ norm of (a limited frequency part of) the fields H and E. These results were obtained for equations posed 
in the whole space (SI = R 3 ) and for constant coefficients e and \x. 

In practice, the various coefficients of the system may not be constant. Typically, the magnetic permeability and 
electric permittivity may depend on the space variable x and have jumps across the boundary of the domain Q. 

Adapting the above mentioned compensated compactness argument, Jochmann established in |22j the existence 
and weak stability of global finite energy solutions for the Maxwell-Landau-Lifschitz system, considering any domain 
SI C R 3 , and variable, possibly discontinuous coefficients (e,fj, e L°°(R 3 )). In the (space) 2-dimensional case, we refer 
to the work of Haddar [20] . 

Concerning the Maxwell-Bloch system, the first author noticed that it shares with the Maxwell-Landau-Lifschitz 
some of its structural properties. This author thus showed in [18] results on existence and uniqueness of global finite 
energy solutions, similar to the ones of Joly, Mctivicr and Rauch, but for some general class of systems including 
the two models above. Again, these results where obtained for equations posed in the whole space and for constant 
coefficients e and \x. 

Here, we continue this study, again for a general class of systems including the Maxwell-Landau-Lifschitz equations 
and the Maxwell-Bloch equations, so as to enlight the similarities and differences between these two models. Adapting 
Jochmann's method, we show the existence and stability of global finite energy solutions, for a given domain SI C K 3 , 
and L°° coefficients. Then, for smooth coefficients, constant out of some compact set, we prove a limit Strichartz 
estimate analogous to the one obtained by Joly, Metivier and Rauch in the constant coefficient case. This allows us 
to show propagation of regularity and uniqueness when initially, curl if and cm\E belong to L 2 (R 3 ). As a corollary 
of a result of Saint- Raymond [38] . we also infer generic uniqueness of the global finite energy solutions. 

Quasi-stationary limits. Next, we turn to the problem of the so-called quasi-stationary limit. Physically, this 
regime appears when the domain SI is small compared to the wavelength. Mathematically, it amounts to some long- 
time asymptotics (replacing in the equations dt by rjdt, for some small parameter rf) with weak nonlinearities (also 
scaled so as to have an amplitude of size n) . 

Jochmann showed in [22] the weak convergence of the corresponding solutions to the Maxwell-Landau-Lifschitz 
system towards the solutions of some reduced system driven by the magnetization, using the weak stability property. 
Starynkevitch extended this result, proving strong and global-in-time convergence in the constant coefficient case in 
|40j . thanks to local energy estimates performed on the explicit fundamental solution of the associated wave equation. 
He also obtained the same result in the case of smooth coefficients, constant out of some compact set, in [5T], thanks 
to dispersive estimates obtained from resolvent estimates on elliptic operators. 

Here, we apply the same methods to our general systems to get weak and strong convergence in the quasi-stationary 
limit. For the latter however, some time integrability assumption is needed to conclude, which is satisfied by the 
Maxwell-Landau-Lifschitz system (since dtM £ L 2 ((0, oo) x SI)), but we do not know if the Maxwell-Bloch system 
enjoys such a property. 

Remark 1. Taking exchange energy into account, one should add to the total magnetic field in ([3]) a term — A" AM. 
The resulting system is then parabolic. We refer to J7|/, J^j, JZ2f, Jiff) /. \ 11$ and for works on the (weak or strong) 
Cauchy problem, and long-time asymptotics. 

2 Main results 

Let us stress that we do not assume that S7 is bounded, for the moment. To deal with both the Maxwell-Landau- 
Lifschitz system ([I])-© and the Maxwell-Bloch system (|4j)-([5]) , we put these two models above into a single class 
of systems consisting in the coupling of the Maxwell equations (with the fields H and E as unknowns) with some 
ODE (corresponding to a third unknown variable). The resulting sytem is symmctrizablc hyperbolic, with semilinear 
nonlincarity, and some structure assumptions arc made, such as affine dependence of the nonlincarity with respect to 
the electromagnetic field, and a priori pointwise estimates on the third unknown variable. One of the key points in 
our study is that the electromagnetic fields decompose into an "irrotational" part, which is directly related to this 
third unknown, and a "divergence free" part, which solves some wave equation. 
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2.1 An abstract setting 

On any finite-dimensional vector space M. N , we denote by u ■ u' the usual scalar product between vectors u and u' , and 
by | • | the associated norm. For all r > 0, B r denotes the (closed) ball centered at 0, with radius r. 
We consider two scalar functions ki(x) and K2(x), which are uniformly positive: 

for i = 1,2, k,G L°°(R 3 ), and 3c >0,k 1 ^ c. (6) 

We denote by H CUI i the space of functions / in L 2 (R 3 ,R 3 ) with curl/ in L 2 (R 3 ,R 3 ). We consider the operator B 
defined by 

B(ui,U2) = (k^ 1 curl U2, — k^ 1 curl ui) for u := (iti,^) G D(B) := i/ cur i x H CUI \. 
This is a skew self-adjoint operator on the Hilbert space L 2 (R 3 ,R 6 ) endowed with the scalar product 

((«l,«2), («j.)«4))ki,k2 := / (niux ■ u[ + k 2 u 2 ■ u' 2 )dx. 



We denote by P(u\,U2) := (PiWi, P2U2) the orthogonal projector on (kerB) 1 - with respect to the weighted scalar 
product above, so that for i = 1, 2, 

ran = {u< e L 2 (R 3 ,R 3 ) | div(«iU*) = 0}, ran (Id -Pi) = {m G L 2 (R 3 ,R 3 ) | curl(uj) = 0}. (7) 

We consider a function F : R 3 x R d x R 6 -> R d , where d <G N, afhne in its third variable, and written 

F(a;,«,u) = F (x,v) +Fi(x,v)u. (8) 

For each j = 0, 1, Fj is measurable with respect to x and continuously differcntiable with respect to v. Furthermore, 



for j = 0, 1, for almost all x € R , Fj(x, 0) = 



(9) 



and Vi? > 0, for almost all x G R 3 , Vi> G + \d v F 3 (x,v)\ < C F (i?). 

Finally, we assume that there exists K ^ such that: 

for almost all x £ R 3 ,V(v,u) G R d x R 6 , F(x,v,u) ■ v ^ K\v\ 2 . (10) 

Remark 2. XTie constant K above may sometimes be taken equal to zero. In this case, Estimate (0) in Theorem\3\is 
improved, since v does not undergo any growth. This is the case for the Maxwell-Landau-Lifschitz model, as well as 
for the Maxwell-Bloch model, when only transverse relaxation is taken into account (Q{p) = —"fPod, for some 7^0, 
and with p oc i the off-diagonal part of p). 

We also consider a function I = (h,l 2 ) G (L°°(R 3 , L(R d , R 3 ))) 2 , where L(R d ,R 3 ) denotes the space of linear 
functions from R d to R 3 . We introduce the following shorthand notation: for any x G R 3 , (k _1 • l)(x) is the mapping 
from R d to R 6 , such that 

for almost all x G R 3 , Vw G R d , (n^ 1 ■ l)(x)v := (n 1 (xy 1 li(x)v, k 2 (x)~ 1 1 2 (x)v). 

Then, for any U := (u, v) in 

L 2 := i 2 (R 3 ,R 6 ) x L 2 (n,R d ), 

the conditions 

div(/tiMi — l\v) = 0, div(K2«2 — fov) = 0, 

may be equivalently written 

(Id - P)(u - (k- 1 ■ l)v) = 0. (11) 
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We look for U £ C([0, oo), L 2 ), with 

v£L^ c ((0,oo),L^(n,R d )), (12) 

solution to 

(d t + B)u = (n- 1 ■ l)F(x,v,u) for i e K 3 , (13) 
d t v = F(x,v,u) for x £ Q, (14) 

and (fTTj) . Here, the solution is understood in the distributional sense, noticing that (fT2]) gives sense to the nonlinear 
term, since the function F(x,v,u) is affinc in u. 

Remark 3. Equations (|ll[) - (|13[) - p4[) reduce to the Maxwell- Landau- Lifschitz system ij7]J-^J when u\ — H , u 2 = E, 
v = M (withd = 3), «i = (i, n 2 = s, h = —fJ>, h = 0, F(x,v,u) = yvA(ui+H a (v)+H ext )—avA(vA(ui+H a (v)+H ext )) 
and to the Maxwell Bloch system when u\ = H, u<z = E, v = p (with d = N 2 ), «i = pi, k 2 = £, h = 0, 

I2 = Tr(r-), F(x, v, u) = —i[A — u 2 ■ T, v] + Q(v). The exterior magnetic field above is usually depending on time. We 
did not consider such time- dependent coefficients in our study, since it would have made notations more intricate; up 
to some integrability assumptions, this extension is straightforward. 

Definition 1. We call U = (u,v) £ C([0, oo),L 2 ) a global finite energy solution to (jli p -Q14 p if (fT2")) holds true and 
U is a solution to pi[) - p4p in the distributional sense. 

Remark 4. Equation (fTTj) has to be seen as a (linear) constraint, which propagates from t = for solutions to 

(HSD-dHD-- 

d t {Id-P)(u- {k- 1 -l)v) =0. (15) 

Indeed, by definition of the projector P, we have (Id — P)B = 0, so that we get (|15p when applying (Id — P) to (|13p . 
using (|14|) (which extends to all x £ M 3 since F(x, 0, u) = 0) and commuting the derivative dt with (Id — P) and k~ 1 - l. 

We therefore have to consider initial data Ui n n satisfying (jlip , and for such constrained initial data, the solutions 
to |J3])-((U]) also satisfy UU as long as they exist. We shall write U mlt := (u init ,v init ) with u init := (u mlt ^,u mlta ). 

Definition 2. Let L div be the set of functions U := (u,v) £ L 2 (R 3 ,R 6 ) x (L 2 (n,R d ) n L°°(n,R d )) satisfying (fTT"j) . 
2.2 Cauchy problems 

Our first result states the existence of global finite energy solutions to (fTTj) - (|T^)l . 

Theorem 3. Assume (|6j) and |8|)- (ll0p . For any Uinit in idiv; there exists U := (u, v) £ C([0, oo), L 2 ), global finite 
energy solution to (jTT|> - (|T4|) wii/i Ui n u as initial data. Moreover, for all T > 0, there is C = C(T,F,l, ||wi„i t ||ioo) such 
that 

(i) for almost all x £ R 3 , for all t^0, \v(t,x)\ < \v in i t (x)\e Kt (with K from ([T0| ); 
(ii) for all t £ [0,T], \\(u,v)(t)\\v < C\\U inU \\v; 

(m) v £ W^((0,oo),L 2 {n,R d )), and for almost every t £ [0,7], \\d t v(t)\\ L 2 (n) C\\U inU \\^. 

Finally, ifUi n it is a bounded set of Ldiv which is compact in L 2 , then for all T > 0, the set U of the above solutions 
with Cauchy data inUi n it is compact in C([0, T], L 2 ). 

To establish this first result, we follow the strategy of Jochmann in [22], which is itself an improvement of the 
method by Joly, Metivier and Rauch in [25 • This is the classical regularization method, in which (global- in-time) 
approximate solutions U n = (u n , v n ) are built first (Section 13. 2p ; the delicate step consists of course in passing to the 
limit n — > oo in the regularization (Section 13. 3p . Pointwise bounds are available for v n , which imply L p bounds for 
(Id — P)u n = (Id — P)(k _1 ■ l)v n , for finite p. The main argument relies on compensated compactness, applied to 
Pu n (Lemma [TBI . 

As a byproduct of the proof of Theorem G2 we also have the following version of stability, where we assume strong 
convergence only for the v part of the initial data. It will be useful below (cf. proof of Theorem [7J when considering 
the weak quasi-stationary limit. 
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Proposition 4. Let {U n ) n m be a sequence in Lf£ c ((0, oo), Lav); bounded in L™ c ((0, oo), L 2 ), im£/i (w n ) rie N bounded 
in W^((0,oo),L 2 (n))nL% c ((0,oo),L°°(ri)) satisfying CEJ), u"|t =0 -> » M zn L 2 (fi) and Bu" = d t D n with (D n ) n 
bounded in L% c {{0, oo), L 2 (R 3 )). Then, up to a subsequence, v n converges to v in L^ c ((0, oo), L p (£l)) for any p ^ 2 
and in Lf^ c ((0,oo),L°°(fl)) weak *, u n converges to u in Lf£ c ((0, oo), L 2 (R 3 )) weak *, and U := (u,v) satisfies (fTT|) . 
dHJ), as we/I as u| t=0 = 

Let us now turn our atttention to smoother solutions. We need to assume more smoothness on the coefficients e 
and fi. Of course, when considering in some physical situation a domain Q with boundaries, the coefficients e and /i 
experiment discontinuity jumps. Since we do not know how to tackle this physical case, we assume from now on that 

n is bounded, and with fC = H, k { - 1 e C^(M. 3 ), i = 1,2. (16) 

In order to get a uniqueness result, we only need to ensure that the "divergence free" part Pu of the fields has the 
H 1 regularity. To this end, once a finite energy solution is given, we make use of the linear system solved by Pu, with 
coefficients depending on the rest of the solution. As in [23], we proceed in two steps: we begin with the propagation 
of regularity, for /i £ (0,1), using Strichartz estimates (Proposition [20]). Applying this result with fj, = 1/2 
provides enough integrability for the coefficients of the above mentioned linear equation to ensure propagation of H 1 
regularity. This implies that u is "almost" L°° , a natural condition to prove uniqueness of the solution. Technically, 
a L°° approximation of u is built thanks to a limit Strichartz estimate for low frequencies (Proposition [2T]) . We also 
need a decoupling assumption, which was introduced in [18] , and is satisfied by the Maxwcll-Landau-Lifschitz system 
as well as by the Maxwcll-Bloch system. 

Theorem 5. In addition to the assumptions of Theorem^ assume (|16p . Let [i £]0, 1], and Umit G L c \- lv with 
curl iti„i ti i £ H^ 1 (M. 3 ), for i = 1,2. Then, the following holds true: 

(a) Any solution U to (|ll|) - (|14p with Uma as initial data given by Theorem^ satisfies curluj £ C([0, oo], H >1 ^ 1 (M 3 )), 
for i = 1,2. 

(b) If [i = 1, assuming moreover that 

there exists j £ {1,2} such that h-jF = and such that F depends only on (x,v,Uj), (17) 
there exists only one solution to (|ll|) - (|14p with Ui n u as initial data as in Theorem^ 

Theorem [5] asserts that the uniqueness property holds for initial data Uinu in Ldw with curliti„i ti i £ L 2 (R 3 ), 
i = 1,2, which are dense in L^w for the topology of L 2 . The following theorem says that the uniqueness property 
even holds generically for the following topologies. Let t s and t w denote respectively the strong and weak topologies 
of L 2 (R 3 ,R 6 ) and let f s denote the strong topology of L 2 (n,M. d ). We consider the product topology r ss (resp. t ws ) 
on L 2 obtained from r s (resp. t w ) and f s . 

Theorem 6. Under the assumptions of Theorem^ [fb]l. for any Ci n u > 0, there exists a Gs dense set L^iv in the set 
{Uinu £ Ldiv | ||uinft||i«(n) ^ Cmit} for the topology t ss and t ws , such that for any U lnlt £ L div , there exists only one 
solution to pi[) - (|14[) with Uinit as initial data, with the same properties as in Theorem^ 

Let us stress that we cannot expect that the problem pT) - p4]) admits smoother solutions than the ones given by 
Theorem [U since, by definition, v is discontinuous across the boundary d£l. However it follows from the general theory 
of discontinuous solutions of hyperbolic semilinear systems [32] [37] that the problem (|TT ]l - (fl~4]l admits piecewise regular 
solutions discontinuous accross 9f2 (let us also refer to the appendix of [H]). Yet the general theory only guaranties 
local-in-time solutions. We do not know if in the particular case of the problem (fTT]) - (fT4"l) global-in-time solutions can 
be obtained. 

2.3 Quasi-stationary limits 

As described in the Introduction, the quasi-stationary regime consists in the limit 77 — >• + for ([TT]) . ([T3]) . (fT4| . where 
dt is replaced with ijdt, and F is replaced with r]F. Equations pT]) and (JTJJ are invariant under this rescaling, whereas 
([TB"]) becomes 

{-qd t + B)u = r^ 1 ■ l)F(x,v,u), for iel 3 . (18) 



6 



For this semi-classical version of (fT4|). it is still true that the constraint (fTTj) is propagated from the initial data. 
Formally, in the limit 77 — > + , w still satisfies (|T4"|). whereas u satifics (|TT|) and _Bu = 0. But for U = (u,v) £ 
C([0, 00), L 2 ), these last two conditions are equivalent to the fact that for all t ^ 0, u(i) is directly determined by v(t), 
and more precisely: 

u= (Id- P)u = (Id- P)(« -1 -l)v. (19) 

Then, f|14|> becomes 

d t v = F(x,v, (Id -P)(k _1 -l)v). (20) 

Using the stability result given by Proposition [4j we have a first result of convergence towards the quasi-stationary 
limit, weakly for u and locally in time for v: 

Theorem 7. Assume (|6j)- (|10p . For any Uinit in Ldw, for any rj <= (0, 1), let U n := (u I? ,u'') be a global finite energy 
solution to ()14p and (|18p with Uinit as initial data. Then, up to a subsequence, v v converges in Lf£ c ((0, 00), L P (Q)) 

for all p ^ 2 and in L^ c ((0, 00), L°°(n)) weak * towards a solution v to (|20|) . with Vinit as initial data; Pu n converges 
to in L^ c ((0,oo),L 2 (O)) weak *, and (Id - P)^ = (Id - P)^ 1 • 0^ converges in Lf£ c ((0, 00), iP(M 3 )) /or aZZ 
p ^ 2 towards u, given by (j 19|) . 

Convergence of the whole sequence t/ 1 ' is ensured as soon as the Cauchy problem associated with the limiting 
equation (|20l) has a unique solution. This is given by the following proposition, which extends |41[ Theorem 3.1] by 
Starynkcvitch. 

Proposition 8. Assume (fTo]) . and letv in i t € L°°(S1). Tftero, there is a unique v G C([0, 00), L 2 (f2))ni^ c ((0, oo),-L°°(0)) 
solution to (|20[) . wiZ/i i^nit as initial data. 

We also prove strong and global-in-time convergence for w, assuming (j!6p again, as well as integrability in time for 
ll^*Hli 2 (fi) an d non-trapping for some wave operator: 

Theorem 9. Under the assumptions of Theorem^ assume moreover (|16p. the non-trapping hypothesis (|1 10|1 and 
that dtv v is bounded (w.r.t. rj) in L 2 ((0,oo) x O). T/ien, Pu v goes to zero in L 2 ((0, 00), L 2 oc (R 3 )). 

Remark 5. In the case of the Maxwell- Landau- Lifschitz system §T§-^§, dtM actually belongs to L 2 ((0,oo) x f2). 
Define the energy £(t) as 

£(t) = l I (e\E\ 2 +fi\H\ 2 )dx+ j fiU(M) + hH ext -M\Adx. 

Differentiating formally this expression with respect to time, we see that the integral of H-cnrlE— E-cutIH = div(EAH) 
vanishes, as well as M ■ dtM (since \M\ is constant). Using the orthogonality relations of the nonlinearity, we get 

d t M ■ H T = a\M A H T \ 2 and \d t M\ 2 = (a 2 + j 2 )\M A H T \ 2 , 

so that estimate (ii) in Theorem^ is improved to 

E ^) + ^T^f j llVM 9 t M ( i ')ll^(fi)^'-^llVM^ex t (i)ll^ W +^ U nM-dtH^dxj dt' = est, 

and the same is true with the quasi- stationary scaling. Assuming for example that H CKt € LfL 2 x and dtH cxt € L\ x , 
we deduce that E is bounded, and dtM belongs to L 2 ((0,oo) x O). 

In the case of the Maxwell-Bloch system, we do not know if such an estimate is available for dtp. 
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3 Existence of global finite energy solutions: proof of Theorem [3] 

3.1 Technical interlude 1 

3.1.1 Intersections and sums of Banach spaces 

We recall some useful properties of the intersection and the sum of Banach spaces. Consider two Banach spaces X\ 
and X 2 that are subsets of a Hausdorff topological vector space X. Then 

x 1 nx 2 :={f e x \ fex u feX 2 } 

(respectively Xx + X 2 := {/ G X | E(f) 0}, where E(f) := {(/i, / 2 ) elixl 2 f x + f 2 = /}) 
is a Banach space endowed with the norm 

H/IUinXa := ||/||xi + \\f\\x 2 
(respectively \\f\\x 1+ x 2 ■= 'm£{\\fi\\ Xl + ||/a||jc a I (/i,/a) € #(/)})• 

If furthermore X x n X 2 is a dense subset of both X\ and X 2 , then (X 1 nX 2 Y = X[ + X' 2 and (Xi + X 2 )' = X[ n 
(cf. Bergh and Lofstrom 0], Lemma 2.3.1 and Theorem 2.7.1). 

3.1.2 Molliflers 

We shall use the following symmetric operators R n : L 2 (M 3 ) — > L 2 (IR 3 ), defined by 

(R n f)(x) := / f(y)w n (x - y)dy for x G M 3 , (21) 

where ai™ G C^°(R 3 ) is a mollifier with supp w n C B(0, 1/(1 + n)) and L 3 iu™ = 1. These operators have the following 
well-known properties: there exists C > such that for all / £ L 2 (M. 3 ), r > 1 and n € N, 

||/ — -R™/||i2(]R3) -> 0, ||-R"/l|i 2 (H3) < C||/|U2(R3), (22) 

\\R n f\\L*(B T )<C\\f\\ LHBr+l) , and ||ii"/|U 2 ( H 3\B r ) < CWfh^B^y (23) 
Moreover for all n £ N, there exists C„ > such that for all / £ L 2 (R 3 ), 

II-R™/I|l~(k 3 ) ^ C n ||/|| L 2(K3). (24) 

3.2 Approximate solutions 

The following lemma claims the existence of global solutions to some regularized problem. 
Lemma 10. For all n G N, there exists U n := (u n ,v n ) G C([0, oo), L 2 ), with 

V n G C([0,oo),L oo (fi,R <i ))nC ,1 ([0,oo), J L 2 (fi,]R d )), (25) 
solution to the regularized problem: 

{d t + B)u n = (k- 1 -/)^ /oriel 3 , (26) 
= F™ /or a; G 0, (27) 

w/iere 

F n (t,x) := F(x,W(t, x),R n u n (t, x)), (28) 
with Uinit as initial data. Moreover, for all n G N, 

(a) For almost all x G K 3 , /or oZi f ^ 0, |w"(<, x)| < \v mlt (x)\e Kt (with K from KTUJ j. 



<S 



(b) ForallT>0, there is C = C(T,F,l, K,\\v in it\\L°°) such that, for all t e [0,T], || (it", v n )(t)\\ L 2 + \\d t v n (t)\\ L 2 {n) <; 

C\\Ui n it\\-L2. 

Proof. The local-in-time solution is constructed via a usual fixed point argument for the mapping A n : C([0, T], L 2 ) — > 
C([0,T],L 2 ), 

= (ex.p(-tB)ui„it + J o cxpdt-t^B)^- 1 ■ l)F n {t\-)dt\ *w + J F> (f , ■)*') , 

where 

:=.F(., *(*,■), '))• 

For T > small enough, A n is shown to be a contraction mapping thanks to properties (|8))- ([T0|) of F, (f24| . and 
because B is a skew self-adjoint operator in the Hilbcrt space L 2 (R 3 ,R 6 ) endowed with the scalar product (■, -) KliK2 . 

Global existence is given by the a priori bounds (a) and (b). The first one follows directly from (fTCTj) and Gronwall's 
lemma. In the same way, taking the L 2 norm of (u n ,v n )(t) = A n (u n ,v n )(t), one gets 

\\(u n ,v n )(t)\\^ «S \\U init \\v+ f (i + ||K- 1 -/|U-)||i ;m (*' ) -)ll^<K / . 

Jo 



One may add to this inequality the one obtained from (|27[) . 

\\d t v n (t,-)\\ L2 ^\\F n (t,-)\\ L2 . 

From (0), j9]), we have 

\\F n (t,-)\\ L 2 ^ C F {\\v m u\\ L ~e Kt )\\(u n ,v n ){t)\\v, 
so that Gronwall's lemma concludes. □ 

3.3 Passing to the limit n — > oo 

Let us stress that Estimate (|24[) is not uniform with respect to n. However, we have: 



Proposition 11. For all T > 0, there is a subsequence of (U n ) ne jq given by Lemma [751 that strongly converges in 
C([0, T], L 2 ) to U := (u,v) G C([0, oo), L 2 ), global finite energy solution to ([Il[) - ()14[) iot£/i t/i„it as initial data, and 
satisfying the estimates §n§, (|m|) o/ 77ieorem[3J 

Proof. First we infer from the bounds (jaj)-((b| in Lemma [TOl that there exists a subsequence, still denoted (it™, v n ), such 
that u" (respectively F n ) tends to u (resp. to Fi im ) in L°°((0, T), L 2 ) weak * (resp. L°°((0, T), L 2 (R 3 )) weak *) and 
v n tends to u in PF^°°((0, T), L 2 (^)) weak * and in L°°((0,T), weak *. This is enough to ensure that 

satisfies (fTT|) . Moreover, Fatou's lemma yields that u and v satisfy (0)-(© of Theorem [3] for almost every t in (0,T) . 

Since the function F is not linear, these weak limits do not suffice to pass to the limit in Equation (|27[) . The 
strategy is to carefully study the nonlinear term F n to prove that the solutions U" of the regularized problems 
(|26)) - ([27|) actually converge (strongly) in L 2 . The key step consists in proving the strong convergence of v n . 

It shall be useful several times to keep in mind that, thanks to the growth conditions (j9|) on F and to the pointwisc 
bound Lemma HU1 (jaj) of the v n , there holds, for all n, m e N, for all (t, x) e [0, T] x M 3 , 

\F?(t,x)\ «S C F (e KT \\v init \\ L ~), (29) 
\F?(t,x)-FT(t,x)\ ^ C F (e KT \\v m . lt \\ L ~W(t,x)-v^(t,x)\, (30) 

where F? = Fi(x,W), i = 0, 1. 

Strong convergence of v n . We perform energy estimates on (|26|) -(|27 |) . Since u may be unbounded, we introduce a 
weight function, which precisely depends on u. More exactly, we choose a positive function po(x) in L°°(R 3 ) n L 2 (R 3 ) 
and define 

p(t,x):=p (x)e- L ti l"^)l ds , (31) 
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with L > C F {e KT \\v lnlt \\L°°)- First, using ([57]) we get 
1 d 



^ t (\\p(v n -v m )\\h m ) (*)= / p 2 K-« m )-(F"-F m )^-L / p^^lK-^l 2 ^. 



Next, decompose f 1 ™ — F m according to (|8]) to get 

In 



~ (IIpK - Olli^n)) (*) = / pV - t- m ) • (iff - K)dx 



(32) 



+ / p (v n — v m ) ■ (F™R n u n — F[ n R m u m )dx — L / p 2 \u\\v n - v m \ 2 dx. 
Jn Jn 

The first term in the r.h.s. of (|3"2"|) can be estimated by C\\p(v n — v m )\\^ 2 ^(t) thanks to (f3T)]) . Now, decompose 
F?R n u n - F[ n R m u m into 

F?R n u n - F™R m u m = F?(R n u n - u) - F[ n {R m u m - u) + (F™ - F™)u. 

The terms produced by the third parenthesis arc estimated thanks to (|30|) . and absorbed by the last term in (|3"2"|) . so 
that 



1 d 

2dt 



(\\p(v n - v m )\\ 2 L2{n} ) (i) < C\\p(v n v m )\\l 2(n) (t) + f p 2 (v n - v m ) ■ F?(R n u n - u)dx 



+ / p 2 (v m - v 11 ) ■ F™{R m u m - u)dx. 
Jn 



Then, decompose R n u n and R m u m according to the orthogonal projector P to get 
1 d 3 

2dT (Jl^-^ll™) (*) < ^IIpK-^II^wW + E^-W + WW' 



(33) 



where 



hi, m , n (t) := / p 2 (v n -v m )-F?R n P(u n -u)dx, 
Jn 



h2,m,n(t) 
h>3,m,n (^) 



p 2 (v n -v m ) ■ F?{R n u-u)da 



The following lemma deals with the term /ii,m,n(i)- 
Lemma 12. There holds 



> 0, 3iV 5 e N, Vn ^ AT 5 , Vm e N, 



C 26. 



(34) 



Proof. First notice that 



J»3 V / 



We first handle the case where x is outside of a large ball. Using the Cauchy-Schwarz inequality, the second 
property of R n in (gU, th e uniform bound in L°°([0, T],L 2 (R 3 )) for u" given in Lemma [TU1 © and the bound (gSJ 
for F™, we get that 



JR 3 \B r 



(j?@*-¥*)F?) ■ P(u n - u)\dxdt \\p(t) 2 \\ L 2 { ^\B r ^)dt 



(35) 



By definition of p there exists r > so that this integral is less than 5. 

It remains to tackle the case where x <G B r . We use the following compactness lemma: 
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Lemma 13 (Jochmann [22], Lemma 3.4). Let (G™)„ g n and (K n ) n eN be bounded sequences in L°°([0, T), L 2 (R 3 , R 6 )) ; 
with K n converging to in L°°([0, T), L 2 (R 3 , R 6 )) weak*. Suppose that (G") ne N is equicontinuous from [0, T] to 
L 2 (R 3 ,R 6 ) and that BK n = d t C n with (C")„ eN bounded in L°°([0, T), L 2 (R 3 , R 6 )). Then for all r > 0, 



sup 



T 

JB T 



G v (t) ■ PK n {t)dxdt 



0. 



(36) 



Let us denote G k > 1 = R k ^p 2 (v k -v l )-Ff) and K n = u n -u. Thanks to ^ and to Lemma[TUl ©, we get that 
(G fe '0fe,/eN and (A'")„ eH are bounded in L°°([0, T), L 2 (R 3 )). Moreover K n tends to zero in L°°([0, T), L 2 (R 3 )) weak *, 
by definition of it. Let us denote T n = J ■ l)F n dt' and T = J (k _1 ■ l)F Vlm dt'. From flU we infer that 

BK" = d t C n , with C n := T n - u n -(J 7 - u). 

The sequence (G n ) n6 N is bounded in L°°([0, T), L 2 (R 3 )). In the same way, equicontinuity is obtained from the bounds 
on dtv n = F n . We therefore apply the lemma observing that, for all m, n G N, 



G m ' n (t) ■ PK n {t)dxdt 



sup 

fc.ZSN 



/ / G k > l {t) ■ PK n {t)dxdt 

JO JB r 



Lemma [T3l therefore ensures that there is N r j G N such that, for n > AV,5 and for all m G N, 



iT^ 2 ^-^)^™) • - u)dxdt 



and Lemma [T2l is proved. 

We now deal with the term /i2,m,n(0- 
Lemma 14. There holds 

^6 > 0, 3N S G N, Vn ^ Vm G N, 



ha, m ,n(t)dt ^6 + C\\p(v n -v m )\\ L 2 \\p(v n -v)\\j 



Proof. The (u n ,v n ) satisfy (| 1 1 [) and so does their weak limit (it, v). Thus 

h 2 ,m.n{t) = - ( p 2 {v n -V m )- F^R n {Id- P)^ 1 -l){~ 1 -v)dx 

Jn 

p 2 (^-^) ■ F{ l R n {Id- P)(k~ 1 • 0(v"-u)da;. 



Then we decompose 



|/l2,m,n(*)|d* < 



(w^ - v^) ■ F? p 2 R n (Id - PXkT 1 • l)(v^ - v) 



pR n (Id - P)p(k~ 1 ■ l)(v n -v) 



dx 



dt 



+ 



pl^-W^) ■ F{ l R n (Id- P) P (k- 1 ■ l)(W*-v)dx 



dt. 



The second integral in the r.h.s. of (|38|) is estimated thanks to Holder's inequality: 



□ 



(37) 



(38) 



p( v n _ B m) . F?R n (Id- P)p(n~ 1 ■ l)(v n - v)dx dt < C||p(i>" - v m )\\ L 2 J\p(v n - ?J)||l? t , 

where C depends only on T 7 F 7 l 7 ||%Mt|U°°- To deal with the first integral in the right-hand side of ([55]). we use the 
following commutation lemma. 
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Lemma 15 (Jochmann [22], Lemma 3.5). Let p belong to L 2 ((0, T), L 2 (R 3 ))nL°°((0, T), L°°(M. 3 )), and let (M")„ eN be 
a bounded sequence in W 1,oo ((0, T), L 2 (R 3 )) n L°°((0, T), L°°(R 3 )) which converges to in L°°((0, T), L 2 (R 3 )) weafc *. 
Then 



T 

2 d"/ 



by 



\\p(tyR n (Id - P)M n - p{t)R n (Id - P)p(t)M n \\ L i {R3)+L 2 {m dt — > 0, (39) 

n— »oo 

||p(t) 2 (W - P)M" - p(t)(Jd - P)/5(t)M™|Ui (K 3 )+i2(K3) ^ — > 0. (40) 
We apply Lemma [TSl (|39|) with M™ = (k _1 • Z)(u n — zJ): the first integral in the right-hand side of (|38| is estimated 



C(T,F,\\v init \ 



f T \\^-^\\ LxnL2 \\ p 2R"(ld - P)M n - pR n (Id - P)pM n \\ Ll+L 2dt, 
Jo 



and thus goes to zero as n goes to infinity, uniformly w.r.t. m. Hence, we get p7[) . □ 

For all S > 0, we also bound ft.3, m , n by <5 for all n ^ N$ and all m € N thanks to (J22J). Finally, summing up with 
(133 and flSTh, we have from (ED: 



V<5 > 0, g N, Vn, m ^ 7V 5 , Vt e [0, T], 

\\p(v n - v m )\\h { n } (t) ^C[6+ \\p(v n - vnWhao^n) + \\p(v n «)||| a((0)T)xn) ) . 

Use Gronwall's Lemma, then let m go to oo, and use Gronwall's Lemma again to deduce: 

V«J > 0, 3N S 6 N, Vn > JV 4l Vt e [0,T], ||p(w" - u)||i»(n)(*) C<5, 

which implies that v n converges towards v strongly in L 2 ((0,T) x ft, p(t 7 x) 2 dtdx). Up to a subsequence, convergence 
then holds almost everywhere, for the measure p 2 dtdx, or dtdx, since p is positive almost everywhere in (0,T) x fl. 
Thanks to the pointwise estimates (jaj) from Lemma ITOl dominated convergence thus ensures that v n converges towards 
v strongly in L 2 ((0, T) x Q, dtdx). Then, equicontinuity of {v} U {i/ l }„ GN in C([0, T],L 2 (Q)) implies (by Ascoli's The- 
orem) the strong convergence of v n in C([0, T], L 2 (VL)). This, together with the uniform bounds on {v} U {u n } n eN and 
with the weak convergence of u n , is enough to pass to the limit in ([2"B]) , ([2~T)l to get (|13p . (HU). 

Strong convergence of u n . Since u" and u satisfy (|26[) and (|13[) respectively, their difference is solution to a hy- 
perbolic equation with source term in L 1 ((0, T), L 2 ' 



(d t + B){u n -u) = (fsT 1 ■ l)(F n - F(x,v, u)). 
The standard energy estimate then gives 

||«" - u\\ L 2(t) [ \\F n - F(x,v,u)\\ L 2(t')dt' 

< C / (||F " - F (x,v)\\ L2 (t') + \\F?R n u n - F 1 {x,v)u\\ L *{t'))dt'. 
Jo 

Thanks to the growth conditions © of F and to the pointwise bound Lemma ITUI (jlj) of the v n , there holds, for 
all n, m e N and (t, x) e [0, T] x R 3 , 

^"(^^-^(x,^,^)! ^ C F {e KT \\v init \\ L oo)\^(t,x)-v(t,x)\. (43) 

In particular this yields that for any t' £ [0, T], \\Fq — Fo(x,v)\\L2(t') goes to zero as n goes to infinity. 
Furthermore, 

\F?R n u n - F x (x,v)u\ < |F"||i?"( u ™ - u)| + I^TIKiT* - Id)u\ + \F? - Fifo v)||u|- 
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Thanks to the L°° bounds on v n (cf. Lemma fTOl (jaj)) and on F[ l (cf. (|29j) ). and to the property (|22|) of the operator 
R n , the first term in the r.h.s. above is bounded by CF(e KT \\vi n i t \\L°°)\u n — u\. In the same way, the second term 
goes to zero in L 2 as n goes to infinity. Finally, up to a subsequence, the third term tends to zero almost everywhere, 
and is bounded by C(F,T, || fi-rtii || ) |w| . By dominated convergence, it thus goes to zero in L 2 . Finally, we get from 
621): 

\\u n -u\\ L *(t)^C(F,T,\\v init \\L°°) I \\u n -u\\ L ,(t')dt' + o(l), 

Jo 

and GronwalPs Lemma shows that u n converges to u in C([0, T],L 2 ). □ 

Thanks to a diagonal extraction process (using times T £ N*), Proposition [TTJ produces U € C([0, oo), L 2 ), solution 
to (fTT ]) -(fT4 l) with Uinit as initial data. Estimates (JTJ) and (jn| are then straightforward. To prove Theorem 02 there 
remains to show its last statement: the stability property. To this end, consider a sequence (Ul nit ) nG ^, bounded in 
£div, and converging to Uinit in L 2 . It generates a (sub)scqucncc of solutions (U n ) n ^, with, from the bounds 
dTTJ) in Theorem [3j ti" converging to u in Lf£ c ((0, oo), L 2 (R 3 )) weak * and v n converging to v in W^^°((0, oo), L 2 (fl)) 
weakly * and in Lf£ c ((0, oo), L°°(Q)) weak *. Then, define the weight p from (j3"Tj) and estimate v n — v m as in (f32j) . 
with R n u n and R m u m replaced with u n and u m , respectively. This leads to the analogue to (j3"3"j) , with no /i3 )TO , n and 
hz, n ,m terms, and no R n in /ii >m , ra and /i2,m,n- Apply Lemma [TBI and Lemma [T5l (HHJ) (instead of (|3T)|) ). to get strong 
C([0, T], L 2 ) convergence of v n towards v (in (|1T|) . the term \\vf nit — v^ it \\ L 2 goes to zero, and contributes to S). Strong 
convergence of the fields u n is then obtained as above , with an initial term ||w"„ it — itimtlli 2 going to zero added to 
the r.h.s. of flU). 

The same process proves Proposition SJ 



4 Propagation of smoothness and uniqueness: proof of Theorem [5] 

It is worth noting that, under the smoothness assumption on e and p in ((T6|) . u G L 2 (IR 3 ,]R 6 ) with Pu 6 i/^M 3 , M 6 ) 
iff u G L 2 (R 3 ,K 6 ) with curlu, G iJ'' _1 (M 3 , R 6 ) for i = 1,2. 

We thus split the proof of Thcorcm[5]in several steps. In Section ETTl we isolate a Cauchy problem for the projection 
Pu of u. This allows some dispersive estimates that we etablish in Section FOl while in Section FL"!?1 Littlewood-Paley 
decompositions are introduced. We consider first the case where p is in (0, 1), then we prove the part (jaj) of Theorem [S] 
in the case p = 1, which concerns the propagation of smoothness, and finally the part (|bj). which concerns uniqueness. 

Remark 6. Let us mention that in the proof of the propagation of H 1 regularity given in fTSjj, the step "p £ (0,1)" 
is missing, and the resulting estimates ( collected here in Lemma [Jity) are claimed without proof. 



4.1 Preliminaries 

Lemma 16. For any solution U := (it, v) to (|11[) - (|14|) with Limit '■— (uinitiVinit) € Ldiv as initial data given by 
Theorem^ the part u := Pu solves for x G M 3 , 

(d t + B)u = P(Au) + Pg, (44) 

U|i= = PUinit, (45) 

where 

A(t,x) := (k" 1 ■ l)Fi(x,v), (46) 

g(t,x) := (k- 1 ■ l)F{x,v, (Id- P)^ 1 ■ l)v), (47) 

= (k- 1 ■ l)F Q (x,v) + (k- 1 ■ /)Fi(x,w)(/d-P)(K" 1 • l)v). (48) 

Proof. First, apply the projector P to the system (fl~3"|) , observing that P commutes with both dt and B. Then, split 

F according to (j8j), split u into u = u + (Id — P)u and finally use the constraint (jTTj) . □ 

The projectors P;, i ~ 1,2, defined on L 2 (M 3 ,IR 3 ), extend to L P (M 3 ,R 3 ) (this result extends the classical one by 
Caldcron and Zygmund [11) on singular integrals, in the spirit of the extension by Judovic [44]): 
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Lemma 17 (Starynkevitch [41] . Lemma 3.13). Under assumption (|16p . the projectors Pi, i = 1,2, extend to F P (R 
and for all po > 1, toere exists C > smc/i that for all p G [po,oo), their norm from F P (R 3 ,R 3 ) into itself are less than 
Cp. 

We deduce estimates for the right-hand side of (|44|) : 

Lemma 18. ^4s in Theorem^ assume ([5|- (|10[) and (| 16[> . Lei t/j„it := (nmit,t>imt) G Pdiv, and let U :~ (u,v) be any 
solution to (| 1 1 1) - (| 14[) with Ui n u '■= {ui n it, Vi n it) G Pdiv as initial data given by Theorem^ The following holds true for 
A and g given by ([4"6 ]) -([l8 |) : 

AeF£ c ((0,oo),F°°(R 3 )), (49) 
A G C([0, oo), F 2 (R 3 )), (50) 
5 t AeL~((0,oo),F 2 (R 3 )), (51) 

.gen Kp<oo C([0,oo),i p (R 3 )), (52) 

d t9 G ni <g<2 L,~ ((0, oo), i p (R 3 )). (53) 

Froo/. For all t, t' ^ 0, there holds 

II^(*)||x,«(rs) «S • ^L-(R3)C F (||t; lmt || L =oe Kt ), (54) 

p(t) - < l^" 1 • l\\L°°(m.3)C F (\\v init \\ L ~e Kt )\\v(t) - v(t')\\ L2(n) , (55) 

\\d t A(t)\\ L 2 (m < ll^- 1 • ZlUc^ajCpdlwinitllL-e^HaewWIIi^n), (56) 



what yields estimates (|49 |) - (|5T |) . 

Since v G C([0, oo), F 2 (ft)) n L% c ((0, oo), £°°(fi)), we have v G n p ^iC([0, oo), F p (f2)) - using the boundedness of ft 
forp < 2, and by interpolation for p > 2. LcmmafTTlthcn yields ([52")). Next, using that |Fj(x,u)| ^ C , F(||i'imt||.L°°e' K ' t )|i;| 
for i = 0,1, we infer from (HD that d t v G ni^ 2 F£ c ((0, oo), F P (R 3 )). Since 

ftflr(t,x) = (k -1 •/){a i ,F (x,iJ) •a t w + F 1 (x,TJ)(/d-P)(K- 1 • i) ■ 9 t u + (^F^x, u) • d t u) • (W - F)(k _1 • /)w}, 

thanks to Lemma [T7l we finally get ([53)) . □ 

Also, for a given u, the "fields part" u is in fact uniquely determined: 

Lemma 19. Let A G F^ c ((0, oo), F°°(R 3 )) and g G £^((0, oo), F 2 (R 3 )). For any u m?;t G F 2 (R 3 ), there exists only 
one solution u G C([0, oo), i 2 (R 3 )) to (|44 p -(|45 j) with Pui n u as initial data. Furthermore, it satisfies u = Fu. 

Proof. Existence is given by Lemma [TBI To prove uniqueness, consider two solutions U! and u 2 in C([0,oo),i 2 
to (|33 ]) -(|4"5 ]) . and T > 0. Then 

Vi > 0, ( Ul - u 2 )(t) = / e' t (*- s ) B FA(u 1 - u 2 )(s)ds, 

Jo 

so that, using (gSJ), for t G [0,T], 

|| (ui - u 2 )(*)|| i 2( R 3) < C(k,T) / || (ui - u 2 )(s)j| i 2 (R 3 ) ds. 



Hence, by Gronwall's Lemma, Ui = u 2 on [0,T], for any T > 0. Thus, there is only one solution u. Finally, in the 
same way, u — Fu simply satisfies: 

Vt^O, (3 t + F)(u-Fu)(f) =0, 
so that u = Fu. □ 
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4.2 Technical interlude 2: Fourier analysis 

We recall the existence of a smooth dyadic partition of unity: there exist two radial bump functions x an d 4> valued 
in the interval [0,1], supported respectively in the ball -B(0,4/3) := {|£| < 4/3} and in the annulus C(3/4, 8/3) := 
{3/4 < \i\ < 8/3}, such that 

V£gR 3 , x(£) + 5>(2-^) = l, V£gR 3 \{0}, ^(2"^) = 1, 

\j - j'\ ^> 2 ^ supp (f>(2~ j -) n supp 4>(2~ j '-) = 0, j ^ 1 => supp x(2~ J '-) n supp (j){2' j -) = 0. 

The Fourier transform T is defined on the space of integrable functions / G i 1 (M 3 ) by (J r /)(0 := L 3 e~ 2t7rx '^ f(x)dx, 
and extended to an automorphism of the space 5'(R 3 ) of tempered distributions, which is the dual of the Schwartz 
space S(R 3 ) of rapidly decreasing functions. 

The so-called dyadic blocks A.,- correspond to the Fourier multipliers Aj := 0(2 _J Z?), that is 

Ajli{x) := 2 3j / h(2 j y)u(x - y)dy for j ^ 0, where h := T~ x dj). 

We also introduce So :— x(D), that is 

Squ(x) := / h(y)u(x — y)dy, where h := J-~ 1 x- 



We will use the inhomogeneous Littlewood-Paley decomposition Id = S-i + 53j'eN A?' which holds in the space of 
tempered distributions 5'(M 3 ), and the homogeneous Littlewood-Paley decomposition Id = X^ez A?> which holds in 
<S^(K 3 ), the space of tempered distributions u such that lim^-oo || Xfc^j Afcii||£oc( K 3-) = 0. 

We now recall the definition of the inhomogeneous (respectively homogeneous) Besov spaces (rcsp. Bp q ) on M 3 
which arc, for A G R (the smoothness index), p, q G [1, +oo] (respectively the integral-exponent and the sum-exponent), 
the spaces of tempered distributions u in 6>'(R 3 ) (rcsp. S^R 3 )) such that 

II/IIb^ 9 (r3) := 1 1 So/ 1 1 1> (»3) + ||(2 :,A ||A : ,-/|| LP(K 3)) J || Z ,( N) (rcsp. \\f\\g^ {M3] := |j(2 jA ||A i /|| £P ( R 3 ) ) : ,-|| i<I ( Z )) 
is finite. These Banach spaces do not depend on the choice of the dyadic partition above (cf. for instance the book 

ED- 

4.3 Dispersion 

Propagation of smoothness or singularities for solutions to hyperbolic Cauchy problems, such as 

Lu := (d t + B)u = f, with u\ t=0 = u iniu (57) 

obeys the laws of geometrical optics. Let us refer here to the survey |19j by Garding for an introduction to the subject. 
The characteristic variety of the operator L is defined as 

Char(i) := {(f, x, t, £) G R x M 3 x C x (M 3 \ {0}) | dot L{t, x, r, £) = 0}, 

where L(t,x,r,^) denotes the (principal) symbol of the operator L, which is the 6x6 matrix 



L(t, x, t, £) = L[x 1 t, £) := Tld + B(x 7 £_), where B(x 7 £) :- 



K^X)- 1 ^ 

K2(x)- l £A- 



For all (x,£) G R 3 x (R 3 \ {0}), the matrix B(x,£) admits three eigenvalues: A±(cc,£) := ±(k 1 k 2 {x)) 1//2 |£| and 0, 
each one with multiplicity 2. The eigenspace associated with the eigenvalue is precisely ran (Id — P). We introduce 

P±(x,0 ■= ^ i L{x,-z,£)- l dz, 

2m J| I -A ± (i,e)|=r 
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where r is chosen small enough for A±(x,£) being the only eigenvalue inside the circle of integration. The matrix 
P±i x ,Q is the spectral projection associated with the eigenvalue A±(x, £), that is the projection onto the kernel of 
L(x, —X±(x, £), £) along its range. These spectral projections are homogeneous of degree w.r.t. £, and the associated 
pseudo-differential operators P± satisfy P+ + P- = Id — P. In addition, Id — P, P+ and P- are orthogonal projectors 
(in the weighted L 2 space introduced in Section |2~T|) commuting with B, and acting on Bcsov spaces. 

Now the point is that considering the Cauchy problem (|57p for solutions u satisfying (Id — P)u = 0, we select the 
branch of the characteristic variety which are curved, what generates dispersion. We shall need the following indices 
pt, ri, q t , si, fi, a and p: 

pi G [2, oo) and 1/n + \jp x = 1/2; (58) 
gi G (1,2] and 1/si + 1/gi = 3/2; (59) 
/i£R, a := n - 1 + 2/pi and p := fj, - 1 + 2/q t . (60) 

Proposition 20. Let pi, r±, q\, s\, fi, a and p be given by (|58I) - (I60[) . Under assumption there is a non- decreasing 
function C : (0, oo) — > (0, oo) such that, for any T > 0, for any initial data Ui n u in such that (Id — P)Ui„u = 0, 
and for any source term f in L Sl ((0, T), B p 2 (^ 3 )) such that (Id — P)f = 0, any (weak) solution u to the Cauchy 
problem (|57[) belongs to L Tl ((0, T), B° 2 (^ 3 )) an d satisfies u = Pu, as well as 

ll-P U llz,'-i((0,T),B£ i 2 (K 3 )) ^ C(T)(\\Puinit\\HK + \\PI\\l s i ((0,T),B^ 2 (R 3 )))- (^) 

In the case p\ — 2 (hence T\ = oo, a — fx), the function u = Pu is even in C([0, T], iJ M (R 3 )). 

The result of Proposition |2"D1 is false for r-y = 2, pi = oo, si = 1, q\ = 2 and n = p = 1, a = 0. However, it is true 
when truncating frequencies. We use, for A > 0, the low frequency cut-off operator S x , which is the Fourier multiplier 
with symbol x\ := x('/ty> where the cut-off function \ G C^°(M. d , [0,1]) takes value 1 when |^| ^ 1/2, and when 
|£| ^ 1. Then, we have: 

Proposition 21. Under assumption there is a non- decreasing function C : (0, oo) — > (0, oo) such that, for all 
A,T > and for any u € C([0, T) , H 1 (R. 3 )) solution to ([57)1 . 

\\S* Pu\\lz((o,t),l^(r3)) < C(T)y/\n(l + XT) (\\d x Pu in it\\ L 2(R3) + \\d x Pf\\ L i^ ^ T ),L 2 (R s ))) ■ (62) 

Estimate (|62|) is proved in the case where the operator P is P := —A in [23] Proposition 6.3. Even in this case 
Estimate (|62|) without the cut-off S x is false [24] , [28] . Proposition [21] extends [23] 's result to (smooth) variable 
coefficients. 

Let us mention that one can deduce from these results some similar estimates for the R^-valued solutions u of 
wave equations of the form: 

(<9 t 2 +A)u = f in M 3 , dZ<u\ t =o = u v for v = Q,l, 

where A := — a(x)A + Y^=i B ) ( x )9 3 + C(x), when a — 1 e C^(M 3 ), a(as) ^ Co > and the Bj and C are in 
C£?(R 3 , A^atxat(R))- These equations stand for the propagation of waves in an inhomogencous isotropic compact 
medium K C M 3 surrounded by vacuum. 

To prove Proposition [201 and Proposition [5T] we use the Lax method, that is an explicit representation of the 
solution which allows to take advantage of oscillations via the method of stationary phase, for each dyadic block, to 
get a pointwisc dispersive estimate. The final step relies on the TT* argument and the summation over the dyadic 
blocks. 

This kind of strategy is now very classical and we refer here to the book [3] by Bahouri, Chemin and Danchin for 
a larger overview of its use and of its consequences. However we did not find Proposition [5T] in the literature so that 
we now detail a little bit its proof. 

Proof of Proposition \2(A Let us first remark that it is sufficient to prove Estimate (|61[) for smooth data (by the usual 
regularization process) and locally in time. More precisely, it suffices to prove that there is a constant C > and 
7\ > such that, for all A > and for all u G C([0, Tx],H lr M?)) solution to ([S7]). 

l|-P U lli»-i((0,Ti),S£ ii2 (]R 3 )) ^ C^PUinitWjjn + \\ Pf\\ L>1 (( ,Ti),B^ 2 (K 3 )) ) ' ( 63 ) 
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Indeed, apply several times Estimate ([63)) on time intervals of the form (kTi , (k + l)Ti), with k an integer ranging from 
zero to the integral part K of T/T\ (plus the interval (KTi,T)): on the right hand side, HP^fcTi)!!^,, is estimated 

by c[\\Pu{{k- l)T x )\\^ + \\Pf\\ L . l ^ k _ l )T u ^)^ i3 {^)))- Summing up gives (JUJ). 

Now, consider the operator S(t) := e~ tB P . It admits a parametrix, and thus is given by a sum S(t) = J+ (t) + 1- (t) 
of operators on the dispersive cigenspaces, which are Fourier Integral Operators: for any smooth function u(x), 

f if *(t.z.O-27H/-f ) 

(I±(t)u)(x) = / eV V Va(t,x,Z)u(y)dtdy, 

and wc drop the subscript ± in the sequel. The phase ^(t, x, £) is real, positively homogeneous of degree one in £, C°° 
for £ 7^ 0, and satisfies the eikonal equation: 

dt*(t,x,t) = ±(k 1 k 2 (x))- 1 / 2 \%(t,x,£)\, (64) 

with 

*| t=0 (x,0=27rx-e (65) 

The amplitude a is in Hormander's class S°, and admits an asymptotic expansion whose successive orders satisfy a 
sequence of linear hyperbolic equations. Such a method was initiated by Lax in his pioneering paper |27j . Because of 
the caustic phenomenon, the lifespan of smooth solutions to Equation (|64[) is limited. However, since "J is homogeneous 
in £ and the set K x S 2 is compact, there exists Ti > such that the solution to (|M1) - (I551 remains smooth on [— T\, T{\. 
Let us mention here that Ludwig |31j succeeded in extending Lax' analysis into a global- in-time result. The arguments 
have been refined thanks to Hormander's theory of Fourier Integral Operators |21j . |17j . But we use the paramctrices 
(and the solution operators I(t)) only locally in time and we refer for their construction to the work of Chazarain [15j . 
Nirenberg and Treves [35] and [M], Kumano-go [33] and Brenner [5]. In particular we refer to the last one for the 
following precious informations about the phase ([§] Lemma 2.1): there exist c, C > such that 

(i) . c|e| < \K\ < on [-Ti.Tj x M 3 x (M 3 \ {0}); 

(ii) . c/d < CId, being real symmetric, on [-T X) Ti] x M 3 x (K 3 \ {0}); 

(iii) . is semi-definite with rank 2 for |£| ^ 0, t 7^ 0; and for |£| = 1, x € K , there is a constant Co > such that 

the moduli of the non-zero eigenvalues of are bounded from below by co\t\; 

(iv) . for x K, the above results are consequence of the exact formula: ^(t, x, £) = 2ttx ■ £ ± 2n(KiK2(x))~ 1 / 2 t\£,\. 

Note that these results imply that the kernel of I(t) is a Lagrangian distribution. 
We use the TT* method for the frequency localized operators 

T 3 {t) := Ajl(t), j e Z. 

The composed operator Tj(i)Tj(t')* is then 

T j (t)T j (t , r = A j I(t-t')A j . 

Here is the TT* result. 

Lemma 22. There exist < c < C such that for all j E Z, u E L 2 (R 3 ), p,r E [l,oo] and f E i r '((0, T),L P ' (M 3 )), 

ll T i(*)^IU'-((o,T),LP(K3)) < ||ui| L 2 (K3) sup (bj(g,g))i, (66) 

= sup 16,(7,5)1, (67) 

L r ((0,T),£*(R3)) se^r.p 



and 



T 3 {t)T 3 {tyf{t')dt' 
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with 



W,9) 



(0,T)x(0,T) 



{T 3 {t)T 3 {t'y f{t') )9 {t)) L 2 (m dtdt' , 



B 3 r p the space of functions g(t,x) in B r ^ p whose Fourier transform is supported in c2 3 ^ |£| ^ C2 J , and B r ^ p the space 
of smooth functions g(t, x) satisfying HffH^T-'^o T) l*>' ^ ^' 

Proof. Begin with 



l|Tj(*) u IU T '((0,T),LJ>(K3)) = Sup 



Tj(t)u{x) ■ g(t,x)dtdx 



'(0,T)xR 3 

Then, use the Plancherel identity plus the properties of the support of the Fourier transform of Tj(t)u to get 



II t jW u IIl''((o,t),lp(k 3 )) = sup 

see* 



/ Tj(t)u-gdtdx 

i(0,T)xE 3 



Using Fubini's principle, transposition and Cauchy-Schwarz inequality, observe that 

i-T 



\\Tj(t)u\\ L 2(( . T ) .l°°(r3)) < IMIl 2 (r 3 ) sup 
Using again Fubini's principle and transposition, get 



T^tfgdt 



(68) 



L 2 (R 3 ) 



= bj{g,g). 



(69) 



L 2 (R 3 ) 



This leads to Moreover, 



T 3 {t)T 3 {trf{t')dt' 



L r ((0,T),LP(R 3 )) 



sup 



(0,T)xR 3 y J0 



TjlfyTjtfYfWdkH) ■ g(t)dtdx 



sup \bj(f,g)\. 



□ 



Now, we need to estimate bj(f,g). We start with a pointwise estimate. 
Lemma 23. There exists C > such that for all j G Z, (t, t') G (0, T) x (0, T) and u G L^R 3 ), 

ll^(*)(^)(*')NU~(R3) ^C2 3 ^(l + 2^|t-* , |)- 1 H il(K 3 ) . 

Proof. Since Aj is a bounded operator on L P (R 3 ) for all p <G [l,oo], it is sufficient to prove that for all t £ (— T,T) 
and u G L^R 3 ), 

||J(t)A i u|| £ o. (B s ) ^C2 3 \l + V\t\)- 1 \\u\\ Ll(m . (70) 
Writing down I(t)AjU, we get, for all x G R 3 : 



(I(t)A jU )(x) = / e W ' ' V Va(t,x,0<P(2- J D)u(y)o%dy 

JR 3 xR 3 



2 3j A e *a**(*.».'?) a ( t>a f > 2> J7 ) v (^)fi(2» J 7) t i» ? , 
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so that 



\(I{t)A jU )(x)\ 2 3j sup 

i/GR 3 



;1 2'(*(t,x,,)-!,.,) a ^ ^ 2 j Tl)<p(T])dTI 



M|z,l(R3). 



To get (|70|) . simply apply the following lemma of stationary phase. 

Lemma 24 (Littman |29j). Let oe a reaZ function C°° such that the rank of its Hessian matrix zs a£ Zeasi p 
and let v(£) be a function supported in a ring. Then there exists M G N and C > (which depends only on a finite 
number of derivatives of ^ , of a lower bound of the maximum of the abolute values of the minors of order p of ^'U, 
on supp v ) such that, for all A G M, 

\\F-\e iA *v)\\ L eo (V)) ^C(l + \A\)- R z- ]T \\D a v\\ LHm . 

| a KM 

To use this lemma, distinguish between short times, for which the eikonal equation (|64[) implies that the phase VP 
admits the expansion 



*(t, x, - 2k(x ■ £ ± «(KiKa(a:))- 1/a |^l) + 0(t 2 ), 



(thus p = 2, A = 2- 7 (kiK2(.t)) 1 l' 2 t is suitable), and subsequent times, for which Estimate (JTTJ) on the phase gives p = 2 
(and A = 2 j ). This yields 



sup 

yGR 3 



< C(l + 2 J 't)" 



(71) 



□ 



Since Tj(t)Tj(t')* is also bounded on L 2 (R 3 ), from the above result and the Riesz-Thorin Interpolation Theorem, 
we infer the following. 

Lemma 25. There exists C > such that for all j G Z, (i,t') G (0,T) x (0,T), p G [2,oo] and u G L P '(K 3 ) 7 

||T,(t)r,(t')^||L P ( R 3) < C2 3 ^ 1 - 2 /rt(l + 2^|i-<'|)-( 1 - 2 ^|| U || iP , (R 3 ) . (72) 

Lemma 26. For any pi,p2 G [2,oo), define r\ and T2 by 1/n + l/pi — 1/2 and l/f2 + l/p2 = 1/2. Then, there exists 
C > sucft i/iai/or a// j G Z, / G L r i ((0, T), M (M 3 )) and 5 G L T '^((0,T), L p ^ (R 3 )), 



Proof. Using (|72|) . we apply Holder's inequality to get 

|6,(/, 5 )| < ^vnll^i + aii.i)-^^!!/!!^,^)^)!!^)!!^^ 



(73) 



L 1 (0,T) 



Using again Holder's inequality, we get 

\b 3 {f,g)\ ^ C2 fl ^i||((i + y|.|)- 3/ri * 11/11^,)) (f) 
Thus, the Hardy-Littlcwood-Sobolcv inequality implies 

^ c ^ / ^ll/ll £ rl( ( o 1 r),£^CK»))'' flr,, i' 1 ((o,r),2^(R»)) 

Moreover, ([69]) yields 



L n(0,T) " 5 "L'''i((0,T),L"i(R3))- 



L 2 (R3) 



= &i(/,/) 1/2 <C2 2 ^ 1 || / ||^ ((oT) ^ (ffi3)) . 



(74) 



(75) 
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Now, we observe that 



bj{f,g)= I (I T 3 {trf{t')dt',T 3 {tYg{t)) L , m dt. 



o jo 

We use the Cauchy-Schwarz and the uniform boundedness of the Tj(t)*, for t in (0,T), in i 2 (M 3 ) to get 

\bj(f,9)\ < bj(f, f) 1/2 \\9\\ LH( o,T),L* m y (76) 

Using now Eq. (|75|) yields 

\b,(f,g)\ C2 2 ^||/|| ir , ((0iT)iil ,, (R3)) || fl || il((0>T)jL2(R 3 )) . (77) 

Interpolating between ([71)) and (1771) . we find (|75|) . □ 
To complete the proof of Proposition [2(3 we need to sum over all frequencies. Apply Duhamel's principle to get 

Pu(t) = e~ tB Pu mlt + f e^-^ B Pf{t')dt'. (78) 
Jo 

The first term in AjPu(t) is thus 

An important fact in order to estimate Ajl(t)ui n it is that the quasi-orthogonality of the dyadic blocks is not destroyed 
by the parametrix I{t). Actually, according to [9l Proposition 1.1 and Lemma 1.4], there exists L £ N* such that for 
all j £ Z, 

A J i(t)= ^imt + Rjit), 

\j-k\*ZL 

where R(t) := X^ez-^iW satisfies 

ll- fl (*) u lli^i((o,r),B; it2 (R3)) ^C\\v\\ 6ll . (79) 
This shows that there is C > such that 

\\AjI(t)Uinit\\LPi(VL s ) ^ C\\Aj I(t)AjUinit\\ Lpi (M 3 ) + || Rj (t)ui n it (K 3 ) , 

so that 

\\I{t) u init\\i J ^i^ t T),i3^ i 2 (R 3 )) 

< C ||||(2 J<T ||A :) /(t)A j M i „ it || L p 1(R 3)) J || i 2( Z )|| LI . i(0/r) + ||||(2 : " T ||i? :) (t)w 4mf || LP i( R 3)) j || i 2( Z )|| Lri(0!T) 

< C\\(2 :, ' J \\AjI(t)AjUi n it\\L r i((0,T),LPi(M 3 )))j\\l 2 (I,) + \\(2 : ' a \\Rj(t)Uinit\\L r i((0,T),LPi(R 3 )))]\\l 2 (Z) 

by Minkowski's inequality 

< C||(2 J '( <H " 2 / ri )|| AjUi n it\\L'(m?))j\\p(Z) + ll(2 JCT ||^i(0 w iriit|l-L r i((0,T),LPi(R 3 )))jlh 2 (Z) 

by J66J and (JJJ) 

^ C||iti n it||^ for some new constant C, using /i = a + 2/ri and ([79")) . 

The estimate for / e [t ' ~ t)B P f(t')dt' follows the same lines, using (|67[) instead of (|66l) . 
Jo 

To conclude the proof of Proposition [20l consider the case pi = 2. As (|6Tj) holds, it is easy to see that u is 
continuous w.r.t. time: using a smooth approximation G L ai ((0, T), if M ) of /, we get Lwfc = and since si ^ 1, 
the usual energy estimates show that (uk)k is a Cauchy sequence in C([0, T], By (|pTj) . as //. tends to / in 

i Sl ((0,T),B^ 2 ), w fe tends to u in C([0, T],H^). □ 
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Proof of Proposition \2l\ Now, in order to prove Proposition [21] it suffices to come back to the proof of Lemma 126] 
taking p\ = oo (and n = 2). Using the standard Young inequality instead of the Hardy-Littlcwood-Sobolcv inequality 
leads to 

< C2 2j ln(l + 2 J T)||/|| L 2 ((0 , T)iL1(R 3 )) ||. g || i2((0iT)iL1(K 3 )) . (80) 
In particular for any g in L 2 ((0, T), L 1 (R 3 )), this yields 

\bj(g,g)\ 1/2 C2VM1 + 2J'T)||s|U a(C o,T),£i(R«)). 

Then we apply Lemma [2"2l with (p, r) = (2, oo) and use (|76|) . where we commute / and to estimate the right hand 
side of dHZJl. We find that for any u £ L 2 (R 3 ), for any / £ L 1 ((0,T),L 2 (M 3 )), 

l|rj(t)u|Ua((o,T),i-(R»)) < C2Vln(l + 2JT)||u|| i 2 (R 3 )) 

|| / r j -(t)r j -(i')*/(* , )^'||L^((o.T).L~(H3)) ^C2Vln(l + 2JT)||/|| L i ((0 ,T),^(H3)). 
Jo 

Then we proceed as in the proof of Proposition [2UJ to sum over the dyadic blocks whose frequencies are below the 
cut-off parameter A, with an extra factor \/ln(l + AT) (and (|75|) holds for n =2, p± = oo). □ 

4.4 Propagation of smoothness: proof of Theorem [5J, case where /i G (0, 1) 

In the case where /i S (0, 1), Theorem [5] is a consequence of the following proposition: 

Proposition 27. Suppose that A and g satisfy the estimates (|49|) - ([53| . that \i £ (0,1) and Uinit £ L 2 (R 3 ,R 6 ) wif/i 
Puimt € if^(R 3 ,R 6 ). Tften i/ie unique solution u e C([0, oo), L 2 (R 3 )) o/ (|34" ]) -([4"5 ]l given 6?/ Lemma [731 belongs to 
C([0,oo),i?^(M 3 )). 

In order to prove Proposition [57] we introduce, for T > 0, the space 

y"(T) := C([0, T], ff"(R 3 )) n C x ([0, T], ff^R 3 )) n L r ((0, T), fl° 2 (R 3 )), (81) 

where p := 2/(l-/i)andr := 2/fx. These indices belong to (2, oo). We shall also use the space (T) := Z?(T)+Z%(T), 
where 

Zf(T) := L 1 ((0,T),^(R 3 ))nC([0,T],L 2 (R 3 )), (82) 
Z%{T) := {/ G C([0, T], L 2 (R 3 )) n L r ((0, T), i?~ 2 (R 3 )) | (83) 
d t f £ L x ((0, T), ^-!(R 3 )) + L s ((0, T), B° 2 (R 3 ))}, (84) 

where q = 2/(2 — fx) and s = 2/(1 + li). These indices belong to (1,2). 

Lemma 28. Suppose that fi, T £ (0, 1), and f £ Z^(T). Then there exists C > such that for any u init in i/ M (R 3 ,R 6 ) 
such that Puinit = Uinit, the unique corresponding solution u £ C([0, oo), L 2 (R 3 )) to Lu = f, with u|t=o = W-init, 
belongs to Y^iT) and 

\\m\\y^ T ) «S C(||ui nit || fff .(]R3) + \\f\\z*(T))- (85) 

Proof. By definition / splits into / = f\ + / 2 , with /, € Z?(T), and u into u = i*o + u i + u 2i where the Ui solve the 
following hyperbolic Cauchy problems: 

Lu = S f with u | t =o = S u init , 

Lui = (Id - 5 )/i with Ui| t =o = {Id - S )u lnlt , 

Lu 2 = {Id - So)fn with M 2 |t=o = 0. 

We already have by energy estimates that the Ui are in C([0, oo), L 2 (R 3 )). Then one gets the estimates of the iij 
in C 1 ([0, T], i/ A '~ 1 (R 3 )) by using the equations. To get the estimates in Y^(T) it therefore only remains to get the 
estimates in L r ((0, T), 5° 2 (R 3 )). 
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For uq this just follows the energy estimate thanks to Bernstein lemma. 

In order to estimate u\ we apply Proposition 1201 with pi = p, n = r, q\ = 2, si = 1, a = 0, p = p. This gives the 
estimate of u x in L r ((0, T), B° 2 (R 3 )) by the right-hand side of flggj). 

In order to estimate U2 we observe that dtU2 satisfies LdtU2 = {Id — So)dtf2, with dtU2\t=o = — <So)/2|t=o- 
By assumption there exists g a G L 1 ((0,T),H^- 1 {R 3 )) and g b G L s ((0, T), B° q 2 (R 3 )) such that <9 t / 2 = .9a + .96- We 
split accordingly <9 t u 2 into <9 t it 2 = u a + u&, where u a solves Lu a = g a , with u a \ t= o = {Id — <So)/2|t=o> and Wf, 
solves Luft = g , with Uf,|t = o = 0. In order to estimate u a (respectively u^) we apply Proposition 1201 with pi = p, 
r\ = r. qi = 2. si = 1, a = — 1, p = p — 1 (respectively with pi = p, n = r, gi = g, si = s, cr = — 1, p = 
and /x — 1 instead of //). This yields the estimate of 9(1*2 in L r ((0, T), i?~ 2 (R 3 )) by ||/2||.z£(T)- As a consequence 
B« 2 = -d t u 2 + {Id - S )f 2 G i r ((0, T), i?~ 2 (R 3 )). Since Pit 2 = u 2 this entails that u 2 G L r ((0, T),B® 2 (R 3 )). 

□ 

The proof of Proposition |2"T1 follows by induction of the following lemma: 

Lemma 29. Suppose that A and g satisfy the estimates (|49p - (|53p . and p G (Oil)- TTien £/iere exists T\ > and 
C > ,suc/i f/ia£ for any Ui n u G L 2 (R 3 ,R 6 ) Pui n it G i? M (R 3 ,R 6 ), i/ie unique corresponding solution u € 

C([0,oo),L 2 (R 3 )) to P4 |) -([45 |l 5 zuen 6y Lemma\M belongs to F^(Ti) and 

ll u l|y"(Ti) < C(||i , «mti||j?M(R3) + |jff|U"(Ti))- (86) 

Proof. In order to prove Lemma [29] we recall the following estimate: 

Lemma 30 (Joly-Metivier-Rauch [23] . Lemma 5.3). There is a constant C , which depends only on ||^4||l«:((o,t),l°o(k 3 ))> 
||9tA|| £ ,«.(( 0l T),i 2 (M 3 )) an d P, such that for all T G (0, 1) and u G Y^T), Aw belongs to Z^{T), and 

\\Au\\ z ^ T ) *S CT^IMIy^t) +C|MIl=°((o,t),l2(R3)). (87) 

Let us warn the reader that there is a small misprint in the right hand side in Lemma 5.3 of |23| . which is corrected 
above. We also have: 

Lemma 31. If g satisfies (|52" ]l -([55 ]) , then for any T > 0, g belongs to Z^{T). 

Proof. Since l<g<2<p<oo there hold continuous embeddings L 9 C B® 2 and L p C B® 2 C B~ 2 , so that, using 
([51 ]) -([53 |1 . we get, for any T > 0, that .g is in Z%(T) C Z^(T). □ 

According to Lemma [30l and Lemma [3~T1 there is a constant C such that for all T £ (0,1), if u G Y^(T) then 
/ := P{Au) + Pg G Z^(T) and 

||/IU"(x) < CT M ^ 2 |M|yM(T) + C|I u IIl°°((o,t),l 2 (e3)). (88) 
Therefore applying Lemma 1281 and choosing Ti small enough we get Lemma 1291 □ 

4.5 Propagation of smoothness: proof of Theorem [5] (jap, case where /i = 1 

We now consider Ui n it in Ldiv with cux\v,i n it,i G L 2 (R 3 ), for i = 1,2, and we consider U solution to (fTT]) - (fl"4"]) with 
Uinit as initial data given by Theorem [3] The idea is to estimate B{Pu) = P(Au) + Pg — dtPu from ([4~4l . 

Lemma 32. Define A and g by (|4l)]) - (|4"5]) . The following holds true: 

a t AGL^ c ((0,oo),L 3 (lR 3 )), (89) 
d tff GL£ c ((0, TO ),L 2 (R 3 )). (90) 

Proof. We apply Theorem[5]in the case p = 1/2. This yields that u belongs to C([0, oo), iJ 1 / 2 (IR 3 )) and thus also to 
I^ c ((0, oo), L 3 (R 3 )). We then infer that d t v G i~ c ((0, oo), L 3 (R 3 )) and then we get the estimates (89) and ([59]). □ 
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Lemma 33. Define A and g by If u G L£ c ((0, oo), H^R 3 )) n W^((0, oo), L 2 (R 3 )) then f := P{Au + g) 

satisfies dtf G I^ c ((0, oo), L 2 (K 3 )), and for any T > 0, there exists C > 0, which only depends on T, A and g, such 
that 

l|5t/||L~((0,T),L 2 (R 3 )) ^ C(!l u llL=°((0,T),ff 1 (R 3 )) + II^ u IIl°°((0,T),L 2 (K 3 )) + 1)- (91) 

Proof. We have dtf — P(dtAu + AdtU + dtg), so that using Holder's inequality, the continuous embedding L 6 (R 3 ) C 
if^R 3 ) and the estimates (|89 ]) -(|90 ]l . we get ((STJ . □ 

Now we observe that <9tu solves for x G R 3 , 

(d f +B)d t u = d t (P(Au) + Pg), (92) 
d t u\ t =0 = -BUinit + P(A| t= oU in it) + Pg\t=o. (93) 

This provides an estimate of d t u in L£ c ((0, oo), £ 2 (R 3 )), and using ([B ]l -([47) ]) . of Bu in L£ c ((0, oo), £ 2 (R 3 )), hence of 
uinL^aO.oo),^^ 



4.6 Uniqueness: proof of Theorem [5] (jbj) 

Let us recall that in this section we assume that there exists j G {1,2} such that h-jF = and such that F depends 
only on (x,v,Uj). Let Umu G Ldiv with curl it j n j ti j G L 2 (R 3 ), for i = 1,2. Let J7 and J7' be two solutions to ([TT]) - (fbI)) . 
given by Theorem [31 both with limit as initial data. The difference SU := U' — U = (Su,5v) between U = (u, v) and 
U' = (u',v') is solution to the following hyperbolic system 

M{5U) := {{d t + B)Su,d t Sv) = ((kT 1 ■ l)SF,SF), where 5F = F(x, v', u'j) - F{x, v, uj). (94) 

Thanks to © we have 

SF = F (x,«',u^) - F (x,v,Uj) + (Fi(x,u') - Fi(x,t;)) • u,- + fj.(x,t/) • (ly - uj)- (95) 

The first and last terms in (|95]) are easily estimated in L 2 (M. 3 ) by C_f (|| i^o || z,=° (K 3 ) e^*) || <5Z7|| i2( K 3) . To deal with the 
second one we construct a L°° approximation of the field uj (analogous to the ones of [23], Lemma 6.2, and [20] . 
Lemma 2.7). 

Lemma 34. There is a non- decreasing function C : (0, oo) — > (0, oo) and for all T > 0, there exists (tin j)\^ e C 
£°°((0,T) x R 3 ) smc/i that for all A ^ e, 

||uj[ >J .(t)||£=.((o,T)xR»)<^(r)lnA ) and || ((Id - - (*)IU»(R») < C(T)/A, /or aii t G [0,T], (96) 
||5 A i^- Wj || jL2((0 , r) , i o C . (R 3 )) < C(T)VhTA, and ||((Jd- 5 A )Pi«j(*)ll£»(R») < ^( r )M, V aH f G [0,T]. (97) 
Let us admit for a while Lemma l34l in order to finish the proof of Theorem[5] Fix T > and consider t G (0, T). In 

the second term on the right-hand side of (|95|) . decompose itj = it| ^. + f (-fd — Pj) u j — u f j) + S X Pj u j + {Id — S x )PjUj. 

We infer from Lemma [34] that 

/ c(r) 

II^IIl 2 (k 3 ) < C F (||fimt||^->(R3)e /f ' r ) ( ||£i/||z,2(R3) + (C(T) In A + \\S x P J u J \\ L ^ iK 3 ) ) ||HU 2 (R 3 ) + 2—7—\\8v\\ L co^ 

(98) 

The energy estimate, together with Gronwall's Lemma, gives 

\\5U(t)\\ L 2 m <: 2C F C(T)jcxp (c F j {l + C{T)\n\+\\S x P jUj {t')\\ L oo^ 3) )dt^j ■ 

Since f*(l + C{T) In A + \\S x P J u :j (t')\\ L o, (¥i 3 ) )dt' < C(T) In A with C(T) we choose T small enough (in order 

to have C(Tq) < 1), and let A go to infinity. This shows that 5U(t) vanishes on [0,To]. Repeat this procedure on 
intervals of size Tq to get the desired conclusion. 



23 



Proof of Lemma \3J\ We define ■ by setting u^j(t,x) := (Id — Pj)uj(t,x) if \(Id — Pj)uj(t,x)\ ^ ClnA, and 
My j(t,x) := otherwise, where the constant C is chosen below (independently of (t,x) £ [0,T] x M 3 ). Therefore, 
for p £ [2, oo), 

wttid-p^-u^mwii = 



(Id—Pj)uj\^C In A 



Kid-p^it^dx^icinxf-mid-p^mi 



(99) 



Now, according to Lcmma fTTl the projection Id—P acts continously in any LP with a norm less than CqP- Furthermore, 
we have 

\\Vi„it\\ P LP ^ \\Vinit\\L°om, 

so that (||i>init ||ip)i<p< 00 is bounded. Thus, using Equation pT|) and the bound from Theorem 131 (jui)) . we infer from 
(1991) that 



IKCJd-.PiK-ufjKt)!! 



(ClnA)P- 2 



(ClnA) A 



2 \2 In 2 



Co ,KT 



choosing p = 2 In A. With C big enough, we obtain (|M| (C(T) = 2C e ifT+1 sup^^ || w in j t ||i ? is suitable). 

We are now concerned with the first inequality in (p?T)) . Coming back to (|44 p -(|45 p . the idea is to use the Strichartz 
estimate (p)2"j) . However, since we are not able to bound d x f in £^ c ((0, oo), L 2 (R 3 )), we cannot apply (|52")) directly. To 
overcome this difficulty we introduce some potential vectors. Since Pu init £ H 1 (R 3 ) (respectively since / = (fi, ,{2) £ 
W^ c °°((0,oc),L 2 (M 3 )) (cf. (Jill) and and P/ = /), for i = 1,2, there exists init := (^,1,^,2) € ff 2 (R 3 ) 

(resp. V := (^1,^2) € ^"((O.ooJ.ff 1 ^ 3 ))) such that 

div(K3_j0»„»t,») = 0, curl(0 lmM ) = KiU in u t i (resp. div(K 3 _ l -0 l ) = 0, curl(^) = 

We consider the operator J3 defined by 

B(<f>i,<t> 2 ) = {H2 1 curl 02,- curl 0i) for ^ : = (01,02) G £>(-B) := D{B) = H cull x # curl . 

The operator B is simply deduced from B by switching hl\ and K2. It therefore shares the same properties and 
estimates. In addition it satisfies the identity: 



n 1 1 curl 






k^ 1 curl 



B = B 



k 1 1 curl 






K7 1 curl 



Let := (0i, 02) be the solution (for x £ R 3 ) of 

(d t + B)(j) = ip, with 0| t= o = cj>init- 
Using the identity (|100[) and Lemma [T9l we obtain 

curl(0j) = KiUi, for i = 1,2. 



(100) 



(101) 



(102) 



Now observe that d t (j) verifies (d t + B)d t (j) = <9tV> , with d t (j)\t=o = "014=0 — B(j) init . Applying the Strichartz estimate 
(l62l) we obtain 



|5 A a t 0|| i 2 ((O T) L oo (R 3 )) sC C(T)vMl + AT)( ||3 x d t 0|t = o||L2(R3) + ||d t d a ,V'l|L 1 ((o,T),L 2 (RS)) 



(103) 



From the definitions of and t/j and from the estimates (|88p and (p?Tj) of the previous sections we get : 

\\d x d t (/)\t=o\\mR3) + \\dtdxi/j\\m(o,T),L2(R*)) < C(T) I] Pu init I] ff i (R 3). (104) 

Using now (|10ip , observing that fe-j and therefore tp3-j vanish because of Assumption [T71 and using (|102p . wc obtain 
the first inequality in (|97p . The second one follows by applying Bernstein lemma. □ 
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5 Generic uniqueness: proof of Theorem [6] 

We apply the following general result of generic uniqueness for evolution equations by Saint-Raymond. 

Theorem 35 (Saint-Raymond [38], Theorem 1). Let £, m -( be a topological space and £ a metric space. Let (S) be an 
evolution equation admitting a solution in £ for any initial data in £i n it- Consider the following hypotheses. 

(HI) For any initial data Ui n u G £i n it, for any {Uf nit ) e tending to Uinit in £init, for any (U E ) e in £ respective solutions 
to (S) with Uf nit as initial data, 

(i) there exists a limit point of {U e ) £ in £ ; 

(ii) any limit point of (U £ ) e in £ is solution to (S) with Uinit as initial data. 

(H2) There exists T>, dense subset of £ , such that for any Uinit in T>, there exists only one solution to (S) in £ with 
Uinit as initial data. 

Under these two hypotheses, there exists a G$ dense £mit of £i n u such that for any Ui n a G £init, there exists only one 
solution to (S) in £ with Uinit as initial data. 

Recall that we denote by r s and t w respectively the strong and weak topologies of L 2 (R 3 ,R 6 ), and by f s the strong 
topology of L 2 (n,M. d ). We consider the product topology t ss (resp. t ws ) on L 2 obtained from r s (resp. t w ) and f s . 
For any Ci n u > 0, consider 

£init := {Uinit G -^div | || v init || L=° (CI) *S Cmit}j 

endowed with the topology t ss (resp. t ws ) inherited from L 2 , and 

£ := {U G C([0,oo),L 2 ), satisfying ([12]) and the estimates (0), flu]), ([m]) of Theorem [3]} , 

endowed with the strong topology (resp. the weak * topology relative to t ws ) of C([0, oo), L 2 ). Hypothesis (E(T]) is a 
direct consequence of the stability property stated in Theorem [3] (resp. Proposition 2]). Now, set 

V := {U init G idiv with curl /;,„,,., G L 2 (R 3 ), for i = 1,2}, 

which is dense in £i n u for the topology t ss inherited from L 2 . Moreover Theorem [S] yields that Hypothesis (1^2) is 
satisfied. We can therefore apply Theorem [3S] what proves Theorem [S] 

6 Quasi-stationary limits: proof of Theorem [TJ, Proposition [8] and The- 
orem [9] 

Proof of Theorem^ We first observe that the bounds (01, (JITJ) given by Theorem [3] for U v are uniform in n G (0, 1). 
Therefore, up to a subsequence, U 71 converges to U := (u,v) in W, '^°((0, oo), L 2 (f2)) weak * and v n converges to v in 
L^ c ((0,oo),L°°(^)) weak *. In addition, there holds BvP = d t D'\ with D r > := -r](u r ' - (k" 1 • l)W). Passing to the 
limit already yields that U satisfies the (linear) equations (fTTj) and Bu = 0. Using Proposition 0] we also get that U 
satisfies (HU, which means that v solves ([2TJ]) . □ 

Proof of Proposition [3 The proof is very similar to the uniqueness proof in Theorem [5] (|b| : it relies on some L°° 
approximation of (Ld-P)(n- 1 ■ l)v. Consider v 1 , v 2 & C([0, oo), L 2 (fl)) n i^ c ((0, oo), L°°(tt)), solutions to ([20]) with 
the same initial data Vmit, and define Sv := v\ — V2- Fix T > 0. From the properties ([8])- (fT0]) of F, we get the pointwise 
estimate 

d t {\6v\ 2 ) ^C F ((\(Id- P)^- 1 -l^l + l^Svl 2 + \(Ld- P)^- 1 -l^vWSvl) on [0,T] x fi, (105) 
for some constant Cp — Cp (||i'mite Js:T ||L°°(o)) • Now, defining for M > 

:= 1 \(Id-P)(K-^l)v 1 \^M(Id - P)^ 1 ■ l)vi, 
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we get from Lemma [34] that there is C(T) > such that 

ViVOl, \\ivf |Uoo ((0>r)xa) ^ M, Wild-P^-l^-wf || L » c( o,D,L»(n)) ^C{T)e~ M ^ T \ 

Integrating (|105[) over f2, using the Cauchy-Schwarz inequality and increasing the constant C (which is still independent 
of M), we obtain 

dt (lN|i a(n) ) < C F ((M + l)||^||l 2(n) +C{T)e- M ' c ^) . 
Then, Gronwall's lemma yields 

Vt g [0,21, \\Sv(t)\\l 2{n) < ^l e CHM + l)T-M/C(T), 

Now, choose T so small that CpMT — M/C(T) < (which is possible, since C is a non-decreasing function of T), 
and let M go to infinity. This shows that Sv vanishes on [0,2]. Repeating the argument on successive time intervals 
yields v\ = v-i- □ 

Proof of Theorem^ For each rj £ (0, 1), consider a solution U v (given by Theorem [3]) to (|TT|) . (|14[) and (fT5)l . Conver- 
gence of v v (and (Id — P)m 7 ') is obtained as in the proof of Theorem [7] above. Now, drop the index r\ for simplicity. 
Then, symmetrizing the system by the change of dependent variables 

~ _ ,1/2 
111 — f^i U'ii 

we get in the distributional sense: 

for % = 1, 2, r/d t Ui + (— l) 3_l .R3_i«3_j = 2 ■ kd t v, 

and therefore, applying dt and combining, 

for i = 1,2, Tj^dfUi + J^iEifii = (-l^iZs-^f Z 3 -;^ + r ? 2 < 1/2 / 4 5 t 2 «, (106) 

where we have set 

z = 1,2, := Kg V 2 curl k~ 1/2 (= #3-i, for the duality in i 2 (R 3 , <£c)). 
System (| 1061) shall be understood as a system of wave equations for the "divergence free" parts KiUi, when 

for i = 1,2, 7T, := k-' P i K i 

Then, tt = (tti,^) is an orthogonal projector in the space L 2 (K 3 , dx) x L 2 (M. 3 ,dx). Furthermore, from the description 
of ran P, ; and ran (1 — Pi) in (0, we deduce that 

for i = 1,2, Riiii — Ri (and ttiR* = R* by transposition). 

Thus, we have finally: 

for i = 1, 2, ifdt-KiUi - Q 1 -K i u l = {-ly^Rs-iK^^h-idtV + rf-Ki^^Ud^v, (107) 

with 

for i = l,2, Q, := -J^Jfc + S 1/2 V (VV_\ div( Kl 1/2 - 

From [HI Lemma 3.10], we know that for i = 1,2, the differential second-order operator (—Qi) is a self-adjoint, 
positive, and elliptic. Thus, with v G L£ c ((0, oo), L°°(Q)) n ^ c °°((0, oo), L 2 (fi)) given, and for given initial data, 

1/2 1/2 1/2 

7TiMi| t=0 = Tj/Cj Minit.i and rj(d t TTiUi)\ t=0 = ryK^ liF{v init ,Ui n it) + {-lyRs-iK^Umu^-i, 
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the solution (7Tiui, 7^2/2) to the linear wave equation system (|107|) is uniquely determined. We recover it via vector 
potentials: defining 

for i = 1,2, tfd 2 ^ - Qi4>% = r)Wi K i 1 ^kdtv, 

U=o = 0, (108) 
Vdt<f>i\ t=0 = ^^|*=o> 
we have 

for f = 1,2, mui = rjdtfr + {-IfR^cp^. (109) 

The problem (jl08p also determines uniquely the vector potentials 4>%- Since Hi4>i also satisfies the problem (|108|) we 
infer that iti4>i = 4>i- Furthermore, UU Lemma 3.10] ensures that for i = 1,2, Qj does not admit as a resonance. 
One then needs to assume the following: 

for i = 1,2, is non-trapping. (110) 

This is enough to apply 

Theorem 36 (Starynkevitch |41j . Theorem 3.2). Let Q be a non-trapping, (L 2 -) self- adjoint, negative, and elliptic 
differential second-order operator, for which is not a resonance. Let s > 1/2, 7 € (—3/2,1/2) and R > 0. Then, 
there exists C > such that: for all (u ,Ui) € i^ +1 (R 3 ) x i^(K 3 ), and / such that {x) s {-Q)~<l 2 f £ L 2 ((0,oo) x R 3 ), 
the solution u to 

d 2 u — Qu — f on (0,oo) x R 3 , with Ml =0 = Uo, dtU\ t=0 = Ui, 

satisfies 

\\(u,d t u)\\ L2{{0 x) ^+i {BR)x £y (BR)) < C (llwoll^+^RS) + HII^RS) + II(^) S (-<3) 7/2 /IIl2((0,oo)xR3)) ■ 

For all /jGl, the space ^(M 3 ) is defined by the norm 

IMIij»(R3) = II(-Q) W2 «IIl2 ( r3). 
We apply the result above to 4>i(rjt,x), whith 7 = and s = 1. This leads to: 

Vi? > 0, 3C fl > 0, ||(0 l ,«9t</>i)llL2 ((OiOCl) ^i (i 3 H)xL 2 (i 3 R)) < C/?(?7 1/2 ||7r ? ;K l 1/2 M mi ^J L 2 (R 3 ) 

+ ??||(a-)7r l K l " 1/2 /j5 t u|| L 2 ((0jOo)ji 2( R 3 )) ). 

The right-hand side is controlled thanks to 

Lemma 37 (Starynkevitch [STj , Lemma 3.11). For aZZ i? > 0, i/iere exists Cr > suc/t £/ia£, if m E L 2 (M. 3 ) and 
supp(m) C Br, then 

for i = 1,2, |7r,m(x)| ^ CR(x}~ 3 \\m\\ 1,2(^3} for a.e.xG R 3 . 
Since 7Ti0i = (fn, by the usual TT* argument, \\Ri4>i\\L 2 (B R ) ^ ll^ill jji( J B R ) ' ano - we deduce from (|109|) : 

fori = 1,2, 71^ = 0(07) m£ 2 ((0,oo),L 2 oc (R 3 )), 

which yields the convergence of Pu v to zero in L 2 ((0, 00), F 2 oc (R 3 )). □ 
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